Entanglement sharing in E (g) e Jahn- Teller model in the presence of a magnetic field 
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We discuss the ground state entanglement of the E ® e Jahn- Teller model in the presence of 
a strong transverse magnetic field as a function of the vibronic coupling strength. A complete 
characterization is given of the phenomenon of entanglement sharing in a system composed by 
a qubit coupled to two bosonic modes. Using the residual /-tangle, we find that three-partite 
entanglement is significantly present in the system in the parameter region near the bifurcation 
point of the corresponding classical model. 
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I. INTRODUCTION 

Tools developed in the realm of quantum information 
theory are increasingly being used to investigate funda- 
mental condensed matter problems In particular, 
many model-systems exhibiting quantum phase transi- 
tions have been explored, and new insights in their be- 
havior has been gained by studying the entanglement , 
the block entropy Q and the fidelity Q. In particular, it 
has been demonstrated in general that, apart from acci- 
dental cancellations, entanglement measures always be- 
come singular near the critical points Q (in the thermo- 
dynamic limit) and exhibit a scaling behavior (for finite 
size systems). 

Moreover, entanglement has been shown to display not 
only the signatures of the critical behavior correspond- 
ing to quantum phase transitions, but also to signal the 
presence of bifurcations in the corresponding scmiclassi- 
cal limit This has been demonstrated, for example, 

in some spin-boson models in the strong coupling regime, 
including the Dicke @ and the Jahn- Teller models [9l[l0|. 
In fact, in the collective Dicke model, the two aspects of 
quantum phase transition and classical bifurcation have 
been shown to be related in the adiabatic limit, in which 
scaling laws have been recently derived for the ground 
state entanglement d, [TJ OjJ . 

Generically, spin-boson models describe the linear cou- 
pling of one [13.fIil.flH or many bosonic modes (typi- 
cally, photons or phonons) with electronic or pseudo-spin 
degrees of freedom, usually represented as two level sys- 
tems (qubits). These models have been used to explore 
environment induced decoherence and have been shown 
to display peculiar properties of entanglement [ItJ • 

In this paper, we concentrate on one model of this 
class, the Jahn-Teller (JT) model pj, [H|, involving an 
electron-nuclei system, in which a doubly degenerate elec- 
tronic state (usually denoted as E) is coupled to a dou- 
bly degenerate nuclear displacement mode (e), with the 
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two bosonic modes coupled to different (orthogonal) spin- 
directions of the qubit. This is one of the most investi- 
gated problems in molecular physics for which a variety of 
interesting quantum properties have been demonstrated, 
despite the fact that the corresponding Hamiltonian is 
not exactly solvable. Particularly relevant from our point 
of view, is Ref. Q , where ground state entanglement has 
been investigated for this model-system in the presence of 
a transverse magnetic field, by making use of an approx- 
imate analytic form of the ground state and of numerical 
diagonalization with a truncated basis. There, by study- 
ing the von Neumann entropy, it has been shown that 
the field forces the coupled system into a maximally en- 
tangled state in the large coupling limit. 

Besides this aspect, the E®t system is interesting from 
many respects. Here we concentrate on its multipartite 
structure. Indeed, the model describes a tripartite sys- 
tem with an Hilbert space structure of the kind 2(g)oo(g)oo 
for which we are able to discuss the sharing properties of 
entanglement in the adiabatic limit. 

In general, quantifying three-partite entanglement is 
an extremely difficult task. For the case of qubits, the 
CKW conjectur e 120 1 , re cently demonstrated by Osborne 
and Verstraete [2ll . |22| |. offers us the powerful instru- 
ments of the monogamy inequality and the residual tan- 
gle, which have been already employed to interpret some 
magnetic behaviors [23] . Related results concerning the 
monogamy have been achieved in Ref. [24J, for the case 
of continuous variables. However, no general method 
has been developed for hybrid systems; that is, those 
including both discrete and continuous variables. These 
systems are extremely interesting for many information 
theoretic applications, including the implementation of 
quantum memories or the possibility of entanglement 
concentration and purification [251 ] . In this respect, we 
think it is interesting to study some relevant case, such 
as the JT model we face in this paper. In a related work, 
Tessier et al. [26| examined the case of two-atom Tavis- 
Cumnmings model, making use of the Osborne formula 
[ItJ to obtain the I-tangle. 

With these motivations, this paper explores the shar- 
ing structure of entanglement of the E ® e JT system in 
the presence of a strong uniform magnetic field, whose 



2 



presence has been shown to give rise to interesting con- 
sequences in connection with the Berry Phase 128]. Our 
approach is based on the adiabatic procedure which has 
been already applied to the case of a qubit strongly cou- 
pled to a single slow resonator [Icl |. 

The paper is organized as follows: In Sec. II we formu- 
late the E ® e model in the presence of a magnetic field 
and discuss its solution in the adiabatic approximation; 
in Sec. Ill various entanglement measures are evaluated, 
for which some analytic approximations are derived in 
Sec. IIVI Finally, Sec. IVl summarizes our main findings. 



where we have introduced the dimensionless parameters 
D = 2A/u and L = 2\/2A/w. 

In the adiabatic assumption of slow bosonic modes, and 
as a first step in the Born-Oppenheimcr procedure, we 
will regard as approximately static and solve the qubit 
dynamics for fixed q. 

More formally, we look for a solution of the bi- 
dimensional Schrodinger equation H\tp) = E\tp) written 
in terms of qubit |x(<? )) an< 4 oscillator ip(q ) functions as 

|V) = / d 2 q\^(q)) = ( d 2 q^(q)\q) ® \ X (q)) (6) 



II. THE E <g> e MODEL AND ITS SOLUTION IN 
THE PRESENCE OF AN EXTERNAL FIELD 



where |x(<f)) are the eigenstates of the "adiabatic" equa- 
tion of the qubit part 



The standard JT model describes a qubit interacting 
with two degenerate harmonic modes (conventionally la- 
belled 9 and e). The model Hamiltonian in the presence 
of an external field is the following 
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where we have chosen unit such that h = c = 1. Here lu is 
the natural frequency of the identical oscillators, A is the 
strength of the magnetic field (taken orthogonal to the 
directions of the couplings) and also represents the qubit 
transition frequency, A is the coupling constant, a = I, 
a x , o~y and a z are the usual Pauli matrix and {qg, q e ) are 
real normal coordinates of the vibrational modes. 

The system is invariant under rotations around the 
magnetic field axis and thus there is a conserved oper- 
ator J z , such that [H, J z ] = 0, and which is given by 



L z a Q + -a z 



(2) 



L z being the z component of the orbital angular momen- 
tum 



qePe - q e po 



(3) 



We will take advantage of this symmetry to employ the 
eigenvalues of J z as labels of the energy eigenstates. 

The ground state of the Hamiltonian will be found 
in the Born-Oppenheimer approximation under the as- 
sumption of a fast qubit, which is easily realized for 
strong external fields (A ^> ui). The whole procedure 
can be followed more plainly by rewriting the Hamilto- 
nian ([T|) in polar coordinates as follows 
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2 g +ql , cj> = arctan {q e /qe)- 



Notice that the qubit dynamics is governed by the effec- 
tive g-parametrized magnetic field 



Q = (Lq cos (j>, Lq sin (f>, D) 



(5) 



e-?|x±(g))=±e(?)|x±(?)) 

which gives the eigenvalues 



Q( q ) = |e| = v/jmy. 

The two eigenstates of Eq. ([7]) are 

|X-(g)) - e-^o(g)|T)-e^ft(9)U) 
|X+(g)) = e-^b{q)\])+e^a{q)\[) 



(7) 



(8) 



(9) 
(10) 



where | |) and | J.) are the ±1 eigenstates of a z , while 



a{q) 




(11) 
(12) 



The eigenvalues can be then considered as distortion of 
the harmonic potential, so that the oscillators are effec- 
tively subject to the adiabatic potentials W± = q 2 ±0(q) 
when the qubit is in |x±)- 

The problem, then, reduces to find the solution of a 
bi-dimensional Schrodinger equation with W as the po- 
tential energy. This is a difficult task, which can be sim- 
plified by exploiting the rotational symmetry. 

Since J z commutes with H and due to the functional 
dependence of the adiabatic qubit eigenstates on the po- 
lar angle </>, we can factorize the oscillator wave function 
in the form 



<p{q,4>) = (27r)- 1 /V J (g)e^ 



(13) 



where j = ±1/2, ±3/2, ... is the eigenvalue of the oper- 
ator J z . 

From Eqs.([5] fTUl) the unitary transformation that diag- 
onalizes the potential energy matrix is obtained as 



U = 



e % ^b(q) e l ?a(q) 
e l ia{q) —e l ib(q) 



(14) 
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The transformed Hamiltonian has the form 

H = U^HU = ~ [(|pf + kf ) a + Q(q)a z + A(q)] 



where 



(15) 



A{q) = U^\p\ 2 U + 2U^pU -p = A cr + A ■ a (16) 
The components of rotated effective field A are 
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and 



A, 



q 2 6(6 + D) 
In the absence of magnetic field (the limit D — > 0). 
1 . 1 



A 



A, 



-L z , A y = A z 







(17) 
(18) 
(19) 

(20) 
(21) 



and the well-known result for the linear E®e Jahn- Teller 
model is recovered, [29j, i.e. 



H = 



2 



d 2 



id_ 

qdq 



- q Co + Lqa z 



(22) 



In the strong coupling limit (L> 1), one can neglect the 
off-diagonal (non-adiabatic) terms in this expression, so 
that the factorization (fT3"| leads to a second-order equa- 
tion for the radial function (pj (q) of two adiabatic poten- 
tial energy surfaces (APES) 



dq 2 



ld_ 

qdq 



q 2 ±Lq 
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^(q)=0 (23) 



where the term j 2 /q 2 plays the role of the centrifugal 
energy. In this case, the ground state is characterized 
by the quantum number j — ±1/2 and is thus doubly 
degenerate. 

The off-diagonal non adiabatic terms can be neglected 
directly in Eq.(|15p under the assumption of a strong 
transverse magnetic field, i.e. D ^> 1. This is the regime 
we will discuss. For comparison, in this limit, the Hamil- 
tonian (fT5l) becomes 
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(24) 
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FIG. 1: Normalized ground state wave function for the oscilla- 
tors in the lower adiabatic potential, for D — 10 and different 
values of a. 



The factorization (|T3|) leads to a different equation for 
the radial function ipj (q) 



dq 2 



d 



l 



i 



ipj(q) = 

with the result that, in the presence of a magnetic field, 
the degeneracy present in the linear JT model is broken. 

When D 3> 1 the motion will remain on the low- 
est Adiabatic Potential Energy Surface (APES) given by 
W- = q 2 — 0(q) and characterized by the quantum num- 
ber j = —1/2 (notice that this implies that the centrifu- 
gal energy equals zero). 

Introducing the dimensionless parameter a = L 2 /2D, 
one can show that for a < 1, the potential W- (q) is just a 
broadened harmonic potential surface with a minimum at 
q = and W_(0) = —D. For a > 1, on the other hand, 
the coupling of the oscillator with the qubit splits the 
lowest APES producing a double-well potential surface 
with (a circle of) minima at 



'D 



« = * = 2 [ a ~a 



1 



with 



W-(go) 



D 
2~ 



(26) 



(27) 



In Fig.|T]), the ground state wave function <^_i/ 2 (g) is 
shown for D — 10 and different values of a. We can see 
that the maximum probability amplitude is always found 
around q$, and that, as a decreases, this moves far and 
far away from the origin. 



III. GROUND STATE ENTANGLEMENT 

The expression of the ground state obtained in the pre- 
vious section enables us to compute the entanglement 
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content of the system. We have three independent sub- 
systems: the qubit, the radial, and the azimuthal degrees 
of freedom in which we have decomposed the two oscil- 
lators (from now on, we indicate these subsystems with 
the labels E, q and 4>, respectively). 

In this section we will evaluate the amount of entan- 
glement for every possible bi-partition and then use the 
monogamy inequality to obtain the residual tangle. First, 
however, we briefly review the formalism employed. 

A. I-Tangle formalism 
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To quantify the entanglement for each of the bi- 
partitions of the model we will make use of the I-tangle 
(301 ] , which for a rank-2 mixed state pab can be explicitly 
evaluated as, [23, 



AB J 



Tr(pABPAB) + 2A^ [1 - Tr(p AB )] (28) 

where par is the result of the action of the universal state 
inverter [301 ] on pab 

Pab = S a ®S b (pab) (29) 

and is the smallest eigenvalue of the M matrix de- 

fined by Osborne [27| which is defined and then evaluated 
for our case in appendix [Al 

The universal inverter Si is defined to map every pure 
state pi = \4>){i)\ into a positive multiple of its orthogo- 
nal projector, i.e. Si(pi) = i>i{I — pi). For an arbitrary 
operator O, it gives 



Si(0)=Vi[Tr(0)I-0] 



(30) 



where z/j is an arbitrary real constant (which we choose 
to be unit). The tensor product in Eg . ([2"9"|) , applied to 
an arbitrary joint density operator pab, is given by 

Sa®S b (pab) = Ia®Ib- pa®Ib-Ia®Pb + Pab (31) 

where pa and pb are the reduced density operators ob- 
tained from pab- Putting everything together, 



MpabPab) = 1 - Mpa) ~ Tr(p|) + Tr(p 



(32) 



For a joint pure state (Tv(p AB ) = Tt:{pab) = 1) the I- 
tangle becomes 

tab = 2 [1 - Tr(p A )] (33) 

where Tr(p 2 A ) = Tr(p%). 

We will employ relations ([2"8|) and ([3"3"|) several times in 
the following. 



B. Ground state density operators 

In our case the ground state density operator takes the 
form 

P = [ d 2 ^V-i/2(9)^ 1/2 (g')l9)(<f'l(lx-(<f))<x-(?')l) 



(34) 



FIG. 2: The dependence of the ground state expectation value 
b z = ((Jz) as a function of the parameter a, for various values 
of D. The solid line corresponds to D — > oo. 



There are six nonequivalent bi-partitions: (i) qubit- 
oscillators E ® (4>q); (ii) angular degree of freedom- 
remainder <f> ® (E q); (iii) radial degree of freedom- 
remainder q®{E4>); (iv) angular degree of freedom-qubit 
(f>®E; (v) radial degree of freedom-qubit q®E; (vi) radial 
degree of freedom-angular degree of freedom (j)® q. 

To start evaluating the various tangles, it is useful to 
re-write the ground state density operator ([3"4"]) as 



P 



where 



l«)l/i)ITXTI(/ilH 
|a)|/i>| DUKMW 



•|&>|/2>U>UK/ 2 |<&| 

|6>|/2>|i><tK/i|<o| 



dqq<p-i/2(<l) a(q) \q) , 
dqqip-i/ 2 {q) Hi) \l) 



(35) 

(36) 
(37) 



are two (non normalized) states of the q-mode, while 
\fi) i — 1,2 are the two relevant (and ortho- normal) 
states of the angular degree of freedom: 



\fl 



2tv 



d<t> 
J2tt 



-itj>\ 



, 1/2 



2 " # 



|0) (38) 



The situation is similar to that described in Ref. [9(: 
the angular degree of freedom is constrained to a two- 
dimensional subspace of its total Hilbert space and our 
tripartite system can be considered as a 2 (g> 2 (g> oo system. 
For the set of states (|36H37p we have: 



1 



where 



b z = - 



wl{<l)^-:dq, 
o 0(<7) 



(39) 
(40) 



is the z-component of the Bloch vector b = (a) and 

f°° La 
(a\b) = (b\a) = J o q<p 2 (q)-^dq = -6 (41) 
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FIG. 3: The equatorial component of the Bloch vector along 
the cf> direction, 6$ = (cos^o-^ +sin0<T H ) shown as a function 
of a, for different values of D. The solid line corresponds to 
D -> oo. 

where = (cos 4>a x + singer,,) is the equatorial compo- 
nent in the <p direction. 

In Figs.® and ©, we show the dependence of the 
ground state expectation values b z and on the dimcn- 
sionless quantity a for various values of the external field 
D (broken lines). The continuous plot describes the case 
of very large field (D — > oo) for which an analytic ex- 
pression is obtained in section IIVI We will see in the 
following that these two parameters completely charac- 
terize the ground state. 

From the plots, one can see that for small interaction 
strengths (that is, small ev's) the external field dominates 
and forces the qubit state along its direction; indeed, b z ~ 
— 1 and &0 ~ 0. On the other hand, for a large enough 
a the qubit is strongly correlated with the angular mode 
4> (loosely speaking, it is 'oriented' along cj>) with a small 
residual polarization along the magnetic field. At a = 1 
a singular behavior is found for very large fields, that is 
analyzed below. 

From eq. (f35|) . the marginal density operators are easily 
obtained. For the partitions <f> ® q and E ® q one has: 

= E < S H S > = M/iXAIH + \b)\M(M{b\ (42) 
p Eq = £ (fi\p\fi) = \a)\ TXT \(a\ + \b)\ \{b\ (43) 

i=l,2 

Tracing over q gives a state for E ® <f> which has a bit 
more involved expression: 

/•OO 

PE<t> = / q(q\p\q)dq 

Jo 

= ^ £ l/i)IT)(TI(/i| + ^|/ 2 )U)U K/al 
+ y (l/i)l T)U KM + l)(T KM) (44) 

As stated above, the reduced density operators are com- 
pletely specified by the three set of states introduced 





1 





D=10 




s 








D=20 








D=30 







6 






G 




D=80 




£■ 

e 















2 







Q 0.5 1 1.5 2 2.5 3 
a 



FIG. 4: The tangle between the qubit and the oscillators 
and between the angular degree of freedom and the rest as 
a function of the interaction strength as measured by a for 
different values of D. 

above for the various sub-systems, and by the param- 
eters b z and b^. 

C. Qubit-oscillators, (^-remainder and q-remainder 
tangles 

In this sub-section, we evaluate the entanglement of 
each one of the three subsystems with the remainder. 
Since the overall state is pure, the procedure is quite 
straightforward. The tangle of the qubit with the two 
oscillators is 

r E{M = 2 [1 - Tr(pi)] =\-b\ (45) 

The tangle between the angular degree of freedom with 
the rest of the system is 

r HEq) - 2 [1 - Tt(pI)] . (46) 

Its expression coincides with T E ^ q ^ since the marginal 
density operator for the 4> degree of freedom has the same 
non-zero entries of the qubit one 

- ^l/iX/il + (47) 

These two tangles are shown in Fig. (0}, where it can 
be seen that the qubit (as well as the <f> sub-system) es- 
sentially factorizes for small interaction strengths. This 
is more and more true for increasing external field and is 
due to the fact that the field itself keeps the spin aligned, 
despite its interaction with the oscillators. For values of a 
larger than 1, the interaction dominates more and more. 
This implies that qubit and angular degree of freedom 
becomes more and more entangled; indeed, the tangles 
saturate to 1 for large enough a's. 

To be more precise, and as better discussed below, the 
ground state contains (for almost every a) essentially bi- 
partite entanglement as these two degrees of freedom cor- 
relate to each other, with very little involvement of the q 
part. 
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FIG. 5: The tangle between between the radial degree of free- 
dom and the remainder as a function of a for different values 
of D. Notice that it is notably different from zero only around 
a = 1. 



To see that this is indeed the case, we start by eval- 
uating the entanglement to which the radial degree of 
freedom participates. The tangle T q {E<s>) is given by 



Tq(E<t>) 



2 [1 - Tr(p 2 q )] ^l-b 2 z ~b 2 4 



(48) 



This function is shown in Fig. (|5|), where one can see 
that the radial degree of freedom is very poorly corre- 
lated with the others. This situation is reminiscent of 
the one obtained when a qubit interacts with a single 
oscillator mode in the presence of a tilted external field 
which gives rise to an "asymmetry" in the adiabatic po- 
tential, see [ioj . The 0-mode, here, plays exactly the 
same role of such an asymmetry. In fact, it destroys the 
correlations between the radial mode and the qubit due 
to the monogamy of entanglement. 

It is noteworthy, however, that the entanglement be- 
tween E and q is more relevant in the region around 
a = 1. Indeed, the maximum of T q (E4>) moves towards 
this point as the held increases and this is exactly the 
point where, in the strict adiabatic limit of very large D, 
the tangle becomes discontinuous. 

We show in the following sections that this is precisely 
the region in which a true three-partite entanglement (as 
measured by the residual tangle) is present. To evaluate 
the three-partite correlations, however, we first need to 
evaluate entanglement for the other possible bi-partitions 
in which one of the three subsystems is traced out. This 
can be done explicitly thanks to the Osborne method 
reviewed above. 



D. Angular degree of freedom-qubit tangle 

After tracing over the radial mode q, the reduced den- 
sity operator for the partition E (g) (f>, eq. (|4"4"|) , can be 
re-written in the form 



1 



PEcf, 



bl + bl 



1 



W){vi\ 



bl 



\V 2 ){V2\ 

(49) 



where 



with 



M = /?il/i)IT) + /3 2 |/ 2 )||) 

\v 2 ) = 7il/i)IT)-7 2 |/ 2 )|l) 



(50) 
(51) 



/3i = 




-1/2 



and 



Ti 



(52) 




-1/2 



> 72 = V 1 - 7i 



(53) 

The vectors \vi),i — 1,2 are the only eigen-kets 
of pe<p with non-zero eigenvalues given by Ti = 

(l± ^Jbl + bl)/2. 

This form (which, by the way, shows that the matrix 
has rank two) is particularly useful to apply the Osborne 
procedure. A straightforward calculation gives the tangle 
in the form 

r** = ^ (l + 2A^ 3 ) + | (l - 2A^) (54) 
where is obtained in appendix [Al 



1 



bl + bl 



(55) 



E. q-(f> and q-E tangles 

The two remaining bi-partitions of the system are those 
consisting of the radial degree of freedom and either the 
angular mode or the qubit. These turn out to have no 
entanglement at all. Indeed, one has: 

1 - b 2 z - b\ 

Tr{p E qPEq) = Tr(p 09/ 5 0(J ) = (56) 



1 - h 2 - h 2 

x (Eq) _ x (4>q) _ 1 °z °4> 
rain rnin 



(57) 



2(1-62) 

Putting everything together in eq. (|2"81) . one has 

TEq = 70g = (58) 

F. Residual tangle 

The amount of entanglement for the various bi- 
partitions that we have evaluated above, do not give by 
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FIG. 6: The I-residual tangle given in eq. ((64}, shown for 
different values of the external magnetic field. 



themselves any indication neither on the sharing prop- 
erties nor on the global, three-partite quantum correla- 
tions. Coffman et al. [20] have explored this problem 
in a system of three qubits and introduced a quantity 
known as the residual tangle, to describe the collective 
entanglement content of a state: 



TABC — TA(BC) ~ T AB — TAC 



(59) 



When subsystems A, B and C are entangled with each 
other, the tangle of A with B plus the tangle of A with C 
cannot exceed the tangle of A with the joint subsystem 
BC. This result has been proved valid for any multipar- 
tite state of qubits [2lj | . 

In the E®e JT model, we cannot simply use the defini- 
tion (|59|) of the residual tangle since our three subsystems 
no longer have equal Hilbert space dimension and sym- 
metry under permutations of the subsystems, which is 
present in eq. (1551 would be lost. 

Tessier et al. [26j have faced a similar problem and 
proposed to generalize the quantity (fS"9"]) by just taking 
the average of the three residual tangles to introduce the 
I-residual tangle which has, by definition, the desired per- 
mutation invariancc: 



where 



TEipq 



r (1) 

1 E<pq 

T (2) 
T E<j>q 

T 0) 
E<pq 



r « +r (2) (3) 
'E<pq T 1 E<Pq ' ' Ecj>q 



T E(<Pq) - T E(j) - TEq 
T (j>(Eq) - T~E(f> — T<f, q 
T q(Ecl>) — TEq - T<p q 



(60) 

(61) 
(62) 
(63) 



In our case, r 



(1) _^(2) 



7^ T E&a> an d one easily obtains 



E<pq 



'E^q 



'E<f,q 



2 ( 

TE^q = ^ 



1 - A 



(E<f>)^ 



(64) 



This quantity is shown in Fig. (|6|), from which the sim- 
ilarity with the plots of Fig. ([5]) can be easily grasped. 
This is due to the fact that the q mode is only involved 



in genuinely three-partite entanglement as it does not 
present any bi-partite quantum correlation neither with 
the qubit nor with the angular mode taken alone. 

Again, we notice that the residual tangle is present 
only within a small region around a = 1. 



IV. ASYMPTOTIC BEHAVIOR OF THE 
ENTANGLEMENT 

In order to obtain an analytic estimation of the physi- 
cal quantities evaluated above and for the various entan- 
glement measures introduced, we would need an expres- 
sion for the ground state wave function ^-1/2(9)- It is 
possible to obtain analytically this function under some 
reasonable approximation for the effective adiabatic po- 
tential. In the following we report three distinct approxi- 
mations, valid in the regimes of i) small coupling, a«l; 
ii) very large coupling a» 1; and Hi) around the cross- 
over value a « 1. 



A. Small coupling regime 

For a <§; 1 the adiabatic potential in the Schrodinger 
equation (|25j) is approximately harmonic, and the main 
effect of the qubit is to re-normalize the value of the 
oscillator frequency by a factor k = y/1 — a. As a result, 
in this regime the adiabatic ground state wave function 
for the oscillator is well approximated by the gaussian 



k 

f-i/2(q) - (2fc) 5 exp<j --q< 



(65) 



By repeating the various steps of the previous section, 
one can obtain approximate expressions for the various 
tangles introduced above, valid to first order in a. For 
example, 



T E(H) 



T 4>{Eq) 



TEqcf, = T q (E<p) 



2a ira 
~D ' TE * ~ 2D 



V 2) D 



(66) 



(67) 



which we checked to be in very good agreement with the 
numerical solution given above, and which describe the 
start-up of entanglement as soon as the interaction is 
switched on. The last equation shows that (to first order 
in a), the radial mode is involved only in three-partite 
entanglement. 



B. Strong coupling regime 

For a 3> 1, the lowest eigenstate should be localized 
at the minimum of the lowest potential surface. There- 
fore, by expanding the potential around this minimum 
[the qo of Eq. (f2"6"]) ] and by retaining up to second order 
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terms, the Schrodinger equation for the lowest sheet can 
be viewed as the equation for a bi-dimensional shifted 
harmonic oscillator. 

Letting q = q — qo, to be the distance from the mini- 
mum, the approximate adiabatic equation for the ground 
state with j = —1/2, becomes 



For a ~ 1 it can be approximated with the quartic ex- 
pression 

v 2 

„2 



- F 2+—[l- — )- F +V - n 2 q 2 + e_y 2 
dq l qa \ Qo J dq 



where vq = ^-{a 2 — 1) is an energy shift, and 



<P-i/a(q) = 
(68) 



= 



a 



1/2 



1 

2^ 



is, again, a renormalization factor for the oscillator fre- 
quency. 

To obtain analytic estimates for large a, we can take as 
an approximate adiabatic ground state for the oscillator 
the wave function 



2\ V 4 

Wo 



exp 



{-f} 



(69) 



In this regime, we have an almost complete quantum cor- 
relation between the qubit and the 4> mode: 



T E(<j>q) = T 4>(Eq) = T E<P 



1 

or 



(70) 



On the other hand, the q mode is almost factorized since 
its wave function is very localized. As a result, the resid- 
ual tangle is very close to zero (the leading contribution 
being of third order in I /a): 



1 



T q{E<j>) 



a 3 D ' 



3a 3 D 



(71) 



C. Critical region 

The coupling value corresponding to a = 1 divides 
an essentially separable regime from an entangled one. 
This point corresponds to a bifurcation in the appropri- 
ate semiclassical analogue and we have shown that 
the region of parameters around a = 1 is the only one 
with non negligible residual tangle. We have also shown 
that, when the magnetic field increases, this cross-over 
becomes more and more sharp until a singular behavior 
is found in the z-magnetization and (as a consequence) 
in the entanglement measures. 

In this section we seek for an analytic description of the 
system in this parameter region and show that a scaling 
behavior is found with respect to D. For this reason we 
call this a critical region. 

Above, we have defined the adiabatic potential as 



W-(q) = q 2 - e(q) = q 2 - ^ + L*q 2 



W-(q) ~ -D + (l-a)q 2 + 



a 
2D 



q 



(72) 



that describes an anharmonic oscillator for a < 1, 
whereas, for a > 1, it is a double- well potential. As in 
the single oscillator case , this implies that a crossover 
between a localized state and a Schrodinger cat-like state 
is obtained. This, in turn, implies a drastic change in the 
behavior of entanglement. 

This approximate potential apparently depends on the 
two independent parameters a and D, but a reduction 
to a single-parametric problem can be obtained with the 
help of Symanzik scaling [3l| . This is done by re-casting 
the Schrodinger equation (always written for j = —1/2, 
see section |lll, into the equivalent form 



dx 2 



1 d 

x dx 



+ (x 2 + x 4 



¥>-i/2(%; C) = e g (C) V-1/2O; C) 



where x = q(a 2 /2Z?) 1 / 6 is a scaled variable. 

remaining scale parameter is, then, £ = {2D j 
a), while the ground-state energy is rewritten as 

f o? \ 1/3 
e ff = e_i /a = -D + I ^=r e g (Q 



The only 



a 2 ) 2/3 (l 



(74) 



It can be shown that all of the qubits and oscillator 
expectation values can be expressed in terms of the di- 
agonal moments: 

f/6 

{q") - I <( v: 1 o ('/)'/</ = | — | {x- 



where 



ciy+ V-i/ 2 ( a; ;C)^ 



(75) 



(76) 



In fact, the parameter £ is very small for a ~ 1 and 
we can obtain analytic approximations for every physical 
quantity we need, by retaining only the first orders of 
their Taylor expansion in 

For example, the two relevant components of the Bloch 
vector of the qubit, taken i) along the external field (b z ), 
and ii) in the equatorial plane along the (f> direction (6^), 
have the approximate expressions 



1 



-V2 



2a 
DP- 

2a 
D^ 



1/3 



1/6 



2a 
£2 



2/3 



(x) 



2a 
D* 



1/2 



(77) 



(78) 



These forms for the components of b can be plugged in 
the general relations for the various tangles obtained in 
section Unl to get 



T E (4,q) ~ T^Eq) ^ 



2/3 



(79) 
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and 



te<i> 



T~Eq<t> = T q(E<j>) 



2/3 



(xf 



2/3 



(80) 



(81) 



All the quantities can be evaluated explicitly once we 
know the various moments of the scaled position x at Q = 
0. These, however, are just constant numerical values, so 
that the physical dependence on D and a can be already 
read from the formula above. In particular, a power- 
law behavior is found, and both the bi-partite and the 
residual tangles become singular as D~ 2 ^ 3 . 

For completeness, we give the numerical values of the 
first moments of the scaled position which are involved 
in the formula above. For a = 1, the problem is reduced 
to the bi-dimcnsional motion in a pure quartic potential 



dx 2 



ld_ 

x dx 



+ x J (p- 1/2 {x; 0) = e_i/ 2 (0)^_x/ 2 (a;; 0) 

(82) 

whose energy and all of the moments can be computed 
numerically. One obtain e_i/2(0) ~ 2.3448, (x) ~ 
0.72737 and (a; 2 ) ~ 0.6515. 

By using these numerical values in Eqs. ([79"]l. (|80[) . and 
(f8Tj) . we obtain that the scaling with D of the various 
tangles at a = 1 is essentially indistinguishable from the 
numerical behaviors for large enough fields (i.e. as long 
as D > 10). 



APPENDIX A: OSBORNE M MATRIX 

The central ingredient required for the computation 
of the /-tangle in Eq. ([28|) is the real symmetric 3x3 
matrix My, derived by Osborne in Ref. [27], for a den- 
sity operator p expressed as a convex combination of its 
eigenvectors: 

P = p\v 1 )(v 1 \ + (l-p)\v 2 )(v 2 \ (Al) 
The independent matrix elements of M are constructed 
in terms of the tensor 

Tijki = Tr(7y7 fe i) 

= Tr(7 i3 )Tr(7 fe; ) - Tr A (Tr B (7y )Tr B (7 fez )) 
- Tr B (Tr A (7 y -)Tr A (7 fei ))+Tr(7„-7 fci ) (A2) 

where 7y = \vi)(vj\. For the partition E®4> one obtains 



Tim 


= 4$$, 






T1112 


= T1121 = 


-2 (A 3 7i - 




T1122 


= T2211 = 


T1221 — T2112 — 1 — 


T1222 


= T2122 = 


-2 {Pill ~ 


Ml) , 


T2222 


= Hil- 







(A3) 



from which, using Eqs. ((52)) and (|53[) . we obtain that the 
only non-zero matrix elements are 



V. SUMMARY 

We have discussed the sharing structure of entangle- 
ment in E eg) e JT model in the presence of a strong 
external field. Using an average residual /-tangle ob- 
tained from the monogamy inequality, we have shown 
that three-partite correlations are important near the 
point in parameter space that corresponds to the bifur- 
cation of the corresponding classical system. This point 
divides a separable from an entangled region, and a singu- 
lar behavior of entanglement is obtained in the strict adi- 
abatic limit. By a detailed analysis performed near this 
point, we have derived a scaling behavior with respect 
to the external magnetic field and identified its "critical" 
exponent. 



M11 



'13 



M33 = 



bl 


bl + 




M31 




t>l- 




bl + 





b z b<p 
bl + bl ' 



(A4) 



The eigenvalues of this M matrix are thus 



A 



{Eft _ 




(A5) 
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